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Abstract 

Antiferroelectric liquid crystals can be considered as a promising al- 
ternative to nematic mixtures in the area of microdisplays. Switching 
behaviour of the molecules has been modelled as two adjacent smectic 
layers. However, some studies made so far reveal that electrooptical 
response can be seriously affected by ion content of the test cells, spe- 
cially in the area of asymmetric cells. Regarding such influence, both 
molecules and ions have been studied in the present work under the 
influence of an external electric field. The mechanisms governing ionic 
behaviour considered so far in the simulations have been diffusion, 
which was dealt through the Nerst-Planck equation, and electrostatic 
interaction, which mixes ionic and antiferroelectric liquid crystal ef- 
fects by means of Gauss equation. In addition, voltage pulses applied 
to the simulated antiferroelectric liquid crystal cell, have been consid- 
ered under different waveforms for possible driving schemes in display 
applications. 

Results show time evolution of ions within the cell, under different 
external applied electric fields and their derealization from alignment 
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layers when high frequency AC pulses are included in the structure of 
the switching waveform. 

1 Introduction. 

Liquid crystals are becoming crucial components in most flat panel 
display applications. Nowadays, lots of nematic mixtures are much 
used to develop a wide range of portable and TV displays with differ- 
ent properties, specially those related to contrast and viewing angle. 
Last advances in research and industrial manufacturing processes are 
overcoming main difficulties with respect to resolution, size, and tem- 
perature range functioning (TJ [2j [3]. Nevertheless, other drawbacks 
have been observed from experimental characterization of market dis- 
plays. This is the case of ghost and still images [I] that can damage 
the display by constituting shadows which are permanent, or have a 
lifetime long enough for being detected by human eye. Such behaviour 
is responsible for a trailing effect clearly observed in video sequences 
5J, and has been related to the slow response time of nematic liquid 
crystals. 

On the other hand, the group of antiferroelectric liquid crystals 
(AFLCs) is being studied and characterized due to their hysteresis, 
experimental fast response time and wide viewing angle. In AFCLs, 
molecules are arranged in layers, called smectic, and can move describ- 
ing two helical motions in the direction normal to the smectic layers, 
with a pitch in the order of 1 — 2{im[6 . Such kind of liquid crystals 
has been often suggested as a good candidate for solving problems re- 
lated to slow response of nematic mixtures, like image trailing in the 
area of microdisplaysO [6l [7j. Hysteresis exhibited by these materi- 
als raises the possibility of developing passive matrix techniques so 
that the cost/quality ratio could be decreased even more than that 
of current TFTs made with nematic mixtures. However, AFLC dy- 
namic properties are experimentally found to be seriously affected by 
changes in the alignment layer. Since these liquid crystals do not have 
a nematic phase within their phase sequence, most commercial poly- 
imides are useless for their alignment 8J. Even those that maintain 
the main characteristics of AFLC response, result in drastic losses of 
contrast, commonly related to a leakage transmission in the absence of 
an applied electric field, called pretransitional effect [9]. 

A summary of the three main possibilities for enhancing contrast 
of antiferroelectric liquid crystal mixtures can be found in [10] . First 
attempts were based on chemical changes that produced a new group 
of antiferroelectric mixtures, called orthoconic liquid crystals 112]. 
and are characterized by their 45° tilt angle for the molecules within 
smectic layers. In surface stabilized structures [HI [14], orthoconic ma- 
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terials behave as isotropic layers, giving rise to the best dark state 
between crossed polarizers [5J [HI US] ■ Apart from changes in the driv- 
ing scheme [T7], other possible solution for improving cell contrast 
can be achieved by producing non-conventional changes in alignment 
layers. An example of this line of research is the sector of asymmetric 
cells, where dissimilar alignment layers are spread over both substrates, 
giving rise to shifted hysteresis cycles [TUBS]. Nevertheless, such hys- 
teresis shift was found to be time-dependent as described in [2D]. The 
main explanation for hysteresis instability lays on ion content of the 
cells, as well as their generation and motion within the alignment lay- 
ers. Ion content mainly comes from synthesis process of liquid crystal 
mixtures [21] or from dissociations impinged by voltage pulses used 
in driving. Ion effects in nematic liquid crystal mixtures have been 
extensively studied by "Liquid Crystals and Photonics" group from 
Ghent University [22j [23] . Their works have shown effects like lateral 
ion transport parallel to alignment layers [HIM], among others. 

The drift-diffusion phenomenon of ions has been also studied in 
other articles, by means of an equivalent three elements circuit [46] . 
Authors relate circuit elements to basic material parameters, specifying 
the validity range of the circuit in each case, with applications in solid 
state electrochemistry systems. 

A model of ion organization in the steady state under the influence 
of an external electric field within a nematic liquid crystal cell, can be 
found in the works of Ionescu, Barbero et al. [26], [27] . They explain 
how switching is enhanced for cathodic polarization while inhibited for 
anodic, in terms of different surface charge densities over LC interfaces 
for both situations. However, nothing is reported about ion evolution 
and molecular movement in dynamic behaviour, that is, when a se- 
quence of different voltage pulses at video frequency is applied to the 
liquid crystal sample under study. 

The present work shows the results of a simulation that consid- 
ers simultaneously ion behaviour and molecular movement under the 
application of certain combinations of voltage pulses. Ions motion 
has been described with the Nerst-Planck mechanism [28], whereas 
molecules have been dealt with a model described in [13j [29] by using 
the Ginzburg-Landau procedure to obtain the motion equations. The 
main objective is to elucidate how changes included in the waveform 
used for driving the liquid crystal cell, can be used for producing ion 
derealization from alignment layers. 

2 Theoretical. 

Main ion sources within a liquid crystal cell have been studied, as well 
as the processes involved in the presence of an applied electric field 
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[301 [3T]. Test cells are supposed to be 1.5/rni thick, and have identical 
alignment layers of nylon 6 over both glass substrates, in a way sim- 
ilar to experimental studies made so far [32, 33J . An antifcrroclcctric 
mixture is considered within such cell in a surface stabilized chevron 
structure [33] [35]. Regarding this arrangement, the displacement of 
cell ion content has been simulated, taking into account phenomenon 
of diffusion under the influence of an external electric field. Simul- 
taneously, molecular switching was considered [T3J [2H] ■ This scheme 
has been applied with different combinations of voltage pulses, in what 
could be a possible waveform for driving the liquid crystal for display 
applications [31)1 1371 155] . 

In the following sections, the main steps for setting this model and 
solving the equations involved are described, as well as the results of 
the simulations performed. 



3 Ionic processes in liquid crystal cells. 

Antiferroelectric liquid crystals contain ions |30j . Main ion sources of 
liquid crystals come from|31j: 

• Liquid crystal synthesis process. 

• Manufacturing protocol of test cells. 

• Molecular dissociation or ion diffusion within alignment layers. 
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Figure 1: Effects on ion content inside a LC cell when an electric field is 
applied: a) Movement; b)Diffusion; c) Formation or annihilation; d) Ad- 
sorption; e) Absorption 



Typical ion concentration in liquid crystal cells oscillates within the 
interval fO 16 — f0 20 ions/m 3 [39]. Therefore, a cell 1/Ltm thick and 1cm 2 
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in area is expected to contain 10 6 — 10 10 ions. Under the influence 
of an external electric field applied to the cell, the ion concentration 
mentioned before evolves being affected by the following combination of 
effects. Some of them are obviously expected regarding arguments from 
classical electromagnetic theory, and are schematically represented in 
the squares of Fig. Q] 

Motion under the electric field: Ion fluxes, J, are represented by: 

J = CqfiE (1) 

Where C is the ion concentration, q the ionic charge , E the electric 
field, and \x the ion mobility. 

Diffusion: This mechanism and motion mentioned before are linked 
by the Nernst-Planck equation [25] : 

dC d 2 n d(nE) 

Where D and fi values are assumed to be constant (i.e. the medium is 
taken as homogeneous). 

Recombination- Generation: Collisions of ions and molecular disso- 
ciations contribute to the temporal evolution of both positive C+ and 
negative C_ ion concentrations, that can be expressed as follows |40j : 

Where, Co is molecular concentration, whereas a and f3 are respectively 
generation and recombination constants. 

Trapping of ions at the alignment layers: Trapped ions form a layer 
over alignment surfaces, so that a screen effect over external electric 
field will be produced. Higher concentrations of free ions (C) within 
the cell are expected to speed up the growth of trapped ions layer. 
Mathematically, this relation could be expressed as follows: 

dt dt t t w 

K is proportional to the trapping probability; C't, C are respectively 
the concentrations of trapped ions and those that have been included 
into the liquid crystal pool. Finally, t t has the meaning of a character- 
istic time of liberation [41] . 

In general, ions generate an electric field against external electric 
field. Thus, they contribute to decrease the external imposed field value 
within the liquid crystal cell. A similar effect happens when AFLC 
polarization is considered. Both effects will be taken into account in 
the model described in the present work. 
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4 One-dimensional ionic model within a 
LC cell. 



A liquid crystal cell behaves as a capacitor, being the antifcrroclcctric 
liquid crystal layer the dielectric material. The glass substrates, cov- 
ered by the semiconductor ITO (Indium-Titanium-Oxide), and nylon 
6 as alignment layer, will act as electrodes for the mentioned capacitor. 
In the next section, a mathematical formalism for ions and molecules 
will be expounded. Main effects under study are diffusion and drift 
caused by the external electric field. 

4.1 Ionic diffusion. 

Ion diffusion can be described in the framework of the Nernst-Planck 
equation [2"B] : 

dC ± _ / D dC ± T c dV\ 
dt dx \ dx dx J 

Being C+ , C_ the concentrations of positive and negative ions, /i the 
ion mobility within the liquid crystal, Co the total ionic concentration, 
D the ionic diffusivity , ze the ionic charge. 
In addition, the following quantities will be defined: 

c= C+ + C_ p= C+_C_ ^_ D 
2Co 2Co ksT 

p is related to the charge concentration. 

C is the total (positive and negative) ionic concentration. 

p is the ionic mobility expressed, by means of Einstein's relation, in 

terms of the temperature T of the cell. 

The Nernst-Planck equations can be rewritten in terms of the new 
quantities p and C as follows: 



dc 






Dze dV 
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(7) 



4.2 Molecular movement. 

In order to describe molecular movement, the model developed by 
Sabater et al. [T5J will be used. This model considers the fol- 
lowing most relevant interactions between molecules: 
Antiferroelectric interaction: 

f c = K V P X ■ P 2 = K, P P 2 S 008(^2 - <Pi ) (8) 



6 



Polarization: 



fpoi = Ps cos 5 I cos v?i — (9) 



Dielectric energy: 

2 

2 



fr>ieX= - \ ^ [ei + Se cos 2 <5 cos 2 ^ + 

Ae(cos sin (5 — sin cos tpi cos (5) 2 ] 



(10) 



av^ 2 

dx 



Elastic energy: 
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2 \cos 2 (5 cos 2 <5 c h 
Total energy: 



(11) 



f = fc + fpol + Jz>ie2 + /jS (12) 

In previous expressions, S is the chevron angle, 5 c i ievron is the 
chevron angle in absence of external electric fields, ifi are the azimuthal 
angles for molecular positions within smectic layers, Ps the material 
polarization, A' , B' , C", I?', L' the liquid crystal elastic constants, and 
£j its dielectric constants. These parameters are shown in Fig. [2] for 
the one-dimensional case. On the other hand, Ae and Se, birefringence 
and biaxiality, respectively, are used to define a new parameter v for 
simplicity within the equations: 

v = Ae sin 2 9 - Se (13) 

Interaction energies (|12[) can be linked by means of Ginzburg- Landau 



df d df 

difi dx ( dtfi 
dx 



(14) 
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Figure 2: Parameters of molecular model in ID. 



4.3 Molecular and ionic interaction: Gauss equa- 
tion. 

Ions and molecules create internal electric fields by the mechanisms 
mentioned so far. Interactions from both sides are considered simulta- 
neously in the equation of Gauss: 

V • D = V (e E + P) = P 2C ze (15) 

Where P as polarization vector is made of Pp, the spontaneous po- 
larization, and Px, a term related to the dielectric susceptibility as 
follows: 

2 

Pp = Ps ^ (— cos tpi cos 6, — sin tpi, cos tpi sin S) 

»=i ^ 
Px = X xx e E 

4.4 Boundary conditions. 

In the process of solving Gauss equation for elucidating ions and molecules 
behaviour in the presence of an external electric field, three effects 
should be taken into account, as boundary conditions: 
Surface interaction, which is considered for molecules in the proximity 



(16) 
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of electrode layers: 



dfo ± dfsurf _ q 



d(fi \ dip 
dx 

, exp [ —asm ] + (1 — c up )exp [ —a cos 



2 J v r \ 2 

I . 2 fi ~ l Pdown\ , . / 2 ^> + <fido 

Cdownexp -asm + (1 - c do wn)exp ' 



fsurf,up,i — 9up 
fsurf,down,i — 9do 

(17) 

Being gi, [3i constants referred to the alignment layers over both up 
and down glass substrates of the test cell [43] . 
Ion confinement within the ce^[28j: 

dc ze dV 

dx ksT^ dx /to\ 
dp _z^ c dV_ _ {j ( ' 

dx ksT dx 

Voltage value in the walls, that should coincide with values externally 
imposed by the driving scheme: 

dV 

V = V±tA s - (19) 
4.5 Simulation parameters. 

Simulations have been done for the AFLC CS-4001 [H]. Parameter 
values for this liquid crystal mixture are shown in Table [I] as well as 
approximations for constants related to surface boundary conditions. 
Some of them were used by Sabater [T31 [221 S3] in his studies for the 
AFLC CS-4000. 



5 Results. 

The model formed by the equations 171 1 21 141 and [T5l and boundary con- 
ditions fTTIfT^l was implemented using a finite differences scheme [IS] . 
Figs. [3H1] show a summary of the results obtained for the simula- 
tions run with different combinations of voltage pulses. In all cases 
ion concentration was fixed to 10 16 ions/m 3 . Information obtained 
has been classified in four groups of results, which include molecules 
and charge evolution data. Regarding external voltage, two waveforms 
were tested: a first bipolar kind and a simple step with a strategy for 
speeding up the relaxation from a switching state to the antiferroelec- 
tric arrangement. Finally, a high frequency AC pulses combination 
was tried to study the effect of fast polarization changes induced by 
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Table 1: Simulation parameters. 



Diffusion 


CS-4001 material constants 


Material constants 


T = 308K 
D = W- u m 2 /s 
q = 1.6 -10- 19 C 
e(Stern layer) = 0.01 
Cq = 10 16 ions/m 3 


A' = 10- n N 
B' = 10- 12 N 

a = \§- 12 n 

D' = W- 12 N 
V = 25000J/m 3 
semi cone = 0.472984rad 
^chevron = 0.418879raa! 
e x = 4.425 • W- u C 2 /Nm 2 
Ae = -4.425 • l$- l2 C 2 /Nm 2 
Se = -8.51 • 10- 13 C 2 /Nm 2 
P s = 79.7 ■ 10" 5 C/m 2 


K p = 10 s 
A = 0.26 
a = 5 
g up = 5 • W- 5 N/m 
9down = 5 • 10- 5 iV/m 

tp up = 2.8274rad 
Vdown = 0.314159rad 
Cup = 0.5 

C-down 0'5 



the external applied voltage in the liquic crystal response. In all cases, 
width and amplitude of pulses used were chosen within the intervals of 
typical switching conditions experimentally observed for the CS-4001 
AFLC mixture in surface stabilized structures [37 . In addition, cited 
waveforms have also shown their utility as possible driving schemes for 
AFLC in microdisplay applications [38]. 

5.1 Simple bipolar waveform. 

The waveform is formed by a pre-switching reset pulse (0 V) for testing 
the antiferroelectric arrangement of molecules (step 1 in Fig. [3]). The 
next pulse is a 25 V selection pulse during 1 ms (steps 2-6 in Fig. [3]). 
Another reset follows the sequence (steps 6-8), and a DC compensation 
for the two last pulses finishes the waveform under study (steps 8-14). 
The data computed from the simulations are grouped as follows: 

• A first group with information about the positions of the molecules. 
It is given through the azimuthal angles, y>\ and (f2, which will 
determine whether the liquid crystal is in one of the ferroelectric 
states (F+, F-) or in the antiferroelectric one (AF). 

• A second group with information about chevron and electric field 
evolution within the cell. 

• A third group with information about the charge concentration, 
p and how it is influenced by changes in the waveform used for 
driving. 
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Thus, it is expected to find AF in the previous reset pulse (step 1 
in Fig. [3]), as well as F+ or F- under the selection pulse influence 
(steps 3-5 in Fig. [3]), that is, the simulation should reproduce such 
results to show an agreement with experimental results reported by 
other authors [37J 13H1 SI] ■ Molecular positions can be observed in Fig. 
|4] for steps 1 and 4 of the waveform proposed in Fig. [3j An analysis 
of the graphs shows that the switching state can be achieved in two 
ways, labeled as a) and b) in the figure. Differences are based on the 
path followed by the molecules to switch from the AF to the F+ state. 
All the simulations carried out show that a) or b) election depends on 
the boundary conditions imposed, as well as on the chevron simulated 
in the middle of the cell. 

For electric charge density, ion migration towards the walls of the 
cell is produced under the switching pulse influence. This charge sep- 
aration is partially maintained during the reset, and can be reversed 
with a switching pulse of opposite polarity. These results are summa- 
rized in Fig. [6] 
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Figure 3: Simple bipolar waveform with reset pulse between selection pulse 
[3-5] and its DC compensation pulse [9-11]. 



5.2 Biasringlet bipolar waveform. 

The waveform is formed by a selection pulse of 35 V during 55/iS (step 
1 in Fig. [7J, followed by a sequence of high frequency AC pulses. Elec- 
trooptical studies made so far [35] show that the period of the AC 
scheme in the bias region (for maintaing the transmission, steps 2-16 
in Fig. [Jj) should be of lower frequency than the one in the section for 
speeding up the relaxation of the AFLC (steps 16-19). A final reset 
at V (step 20) completes the blanking sequence for leading the liq- 
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Figure 4: Plots of the molecular azimuthal angles in the non-switched (step 
1, left), and the switched (step 4, right) configurations. In the last case, the 
ferroelectric state was reached in two ways (a and b). 




AF state F state AF state F state 

Switching a) Switching b) 



Figure 5: Switching process from the antiferroelectric to one of the ferro- 
electric states. Cases a) and b) are the same, although mathematically they 
can be considered as different situations. 



uid crystal to the antiferroelectric state. Regarding ion behaviour, fast 
changes of polarity are expected to produce an effective derealization 
from alignment layers. Fig. [S] shows an evolution of the net charge 
concentration obtained from our simulations, in the main steps of the 
waveform depicted in Fig. [7J It has been theoretically proved that 
fast changes of polarity contribute to charge attenuation and ion der- 
ealization from alignment layers. These results were only found when 
molecular movement was considered, so can be regarded as a result 
of the interaction between ions and molecules in switching conditions. 
An analysis of the plots for azimuts ipi and ip% reveals that the cell 
becomes partitioned in multiple regions with different switching be- 
haviour, as is shown in Fig. [5] This kind of response to external pulses 
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Figure 6: Charge density within the cell in different steps of the bipolar 
waveform applied. Information of the time period is also shown in each 
graphic: a) Step 3 (1100 lis); b) Step 6 (2100 /is), step 8, (3000 /as), and 
step 9 (3100 /is). This plot shows progressive charge derealization from 
alignment surfaces, c) Step 10 (3500 fxs). 



can create a variety of grain boundaries which contributes to ion dis- 
persion specially observed in the plots b), e) and h) of Fig. [H In fact, 
the mentioned derealization seems to be promoted by an increase in 
the frequency used for changing the polarity of external voltage pulses, 
accounting for an attenuation of around a 1% of the total charge in- 
volved. Molecular relaxation is also faster than that achieved with a 
simple bipolar waveform, as has been reported by electrooptical studies 
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Figure 7: Biasringlet bipolar waveform, characterized by secuences of pulses 
of different polarity in both bias and relaxation sections. 
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6 Conclusions. 



A model of a surface stabilized antiferroelectric liquid crystal cell has 
been developed that takes into account molecular and ionic motions in 
the presence of an external electric field. Simulations were performed 
for different combinations of voltage pulses, in what have been con- 
sidered valid waveforms for driving these cells in electrooptical studies 
made so far p~7l [18l QjH [20] . The results first reproduce switching con- 
ditions when external voltage achieves threshold and saturation levels. 
Ion motion creates an internal field opposed to the external one. How- 
ever, fast changes of polarity in the external field produce ionic der- 
ealization from alignment layers as well as weak net charge attenuation 
for high frequency AC pulses. Fast changes of external electric field 
can produce unstable molecular electric fields that interact with ions. 
These local fields produce ionic concentration bursts "far away" from 
the alignment layers. As this effect does not appear in simulations 
developed in the absence of molecular movement, ionic derealization 
has been related to the dynamical interplay between ions and molecules 
under high frequency AC pulses. As a result, driving schemes based 
on fast changes of polarity can be useful in the task of controlling side 
effects of ion concentration within AFLC cells, as possible delays of 
the LC response against external voltage pulses, or cell aging caused 
by driving. 

Further studies with other waveforms and different LC constants 
are under study. 
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Figure 8: Evolution of the net charge density within the cell in strategy 
steps of the Biasringlet waveform. The step is marked in the upper right 
corner of each plot, followed by the time elapsed so far in microseconds. 
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